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Abstract: In supersymmetric models with a long-lived stau being the lightest Standard 
Model superpartner, the stau abundance during primordial nucleosynthesis is tightly con- 
strained. Considering the complete set of stau annihilation channels in the minimal super- 
symmetric Standard Model (MSSM) with real parameters for scenarios in which sparticle 
coannihilations are negligible, we calculate the decoupling of the lighter stau from the 
primordial plasma and identify processes which are capable to deplete the resulting stau 
abundance significantly. We find particularly efficient stau annihilation at the resonance of 
the heavy CP-even Higgs boson and for a lighter stau with a sizeable left-right mixing due 
to enhanced stau-Higgs couplings. Even within the constrained MSSM, we encounter both 
effects leading to exceptionally small values of the resulting stau abundance. Prospects 
for collider phenomenology are discussed and possible implications of our findings are ad- 
dressed with emphasis on gravitino dark matter scenarios. 
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1. Introduction 

The appearance of the lighter stau t\ as the lightest Standard Model superpartner — or 
lightest ordinary superpartner (LOSP) — is a commonplace occurrence even in supersym- 
metric (SUSY) models with restrictive assumptions on the SUSY breaking sector such as 
the constrained minimal super symmetric Standard Model (CMSSM). If the lightest su- 
persymmetric particle (LSP) is assumed to be the LOSP, this parameter region is not 
considered because of severe upper limits on the abundance of massive stable charged par- 
ticles [1]. However, for example, in axino/gravitino LSP scenarios [2, 3, 4, 5, 6] and in 
scenarios with broken R parity [7, 8, 9, 10, ll], 1 the f LOSP becomes unstable and thereby 
a viable option. Indeed, supersymmetric models with a long-lived t\ LOSP are particu- 
larly promising for collider phenomenology [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]: 
Since the t\ LOSP could escape the collider detector as a quasi-stable muon-like particle, it 
can be associated with signatures that are very different from the excess in missing energy 
expected in neutralino LSP scenarios. 

1 In this work we assume that R-parity is conserved. 
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In the early Universe the negatively charged LOSP ri's and the associated positively 
charged anti-staus rj^'s were in thermal equilibrium for temperatures of T > m^/20 > Tf. 
At Tf, the annihilation rate of the (by then) non-relativistic ri's becomes smaller than the 
Hubble rate so that they decouple from the thermal plasma. Thus, for T < Tf, their yield 
Y~ = (n^ + n^*)/s is given approximately by 1~ « Y~ q (Tf), where n~ q ^ = + 
is the (equilibrium) number density of both t\ and and s = 2ir 2 g*s T 3 /45 the entropy 
density with g*s effective degrees of freedom. This thermal relic abundance is subject 
to cosmological constraints in SUSY scenarios with a long-lived t\ LOSP: 

• In axino/gravitino LSP scenarios, governs the non-thermally produced (NTP) 
relic density of axino/gravitino dark matter that originates from T\ decays [3, 5] 
Q^Jj^h 2 = m~^Yf s(To)h 2 / p c where m^/Q denotes the axino/gravitino LSP mass 
and p c /[s{To)h 2 \ = 3.6 x 1CU 9 GeV [1]. Thus, the dark matter density fidm which 
limits U^/^? from above implies an upper limit on Y~ for a given m~,p^. This limit can 

a/G a / ( - r 

become particularly restrictive in the case of additional sizeable contributions to O^m 
such as the ones from thermal axino/gravitino production [25, 26, 27, 28]. For 

example, for = 50 GeV and fij^ = 0.99 S7 dm (0.9$7 dm ), one finds Y~ < 10 -13 

(10- 12 ); cf. Fig. 13 of Ref. [22]. 

• For ri decays during/after big bang nucleosynthesis (BBN), the Standard Model 
particles emitted in addition to the axino/gravitino LSP can affect the abundances 
of the primordial light elements. This leads to upper limits on £ e m/had = e cm/had 
that depend on the stau lifetime [29, 30, 31]. Here e cm /had denotes the (average) 
electromagnetic/hadronic energy emitted in a single f\ decay, which can be calculated 
with particle physics methods for a given model. Accordingly, the BBN constraints on 
Cem/had can be translated into upper limits on 1~; cf. Fig. 12 of Ref. [22] (and Figs. 14 
and 15 of Ref. [32]) for associated limits in gravitino LSP scenarios, which can be 
as restrictive as 1~ < 10 -14 (10 -15 ). 

• The mere presence of the negatively charged ri's at cosmic times of t > 5 x 10 3 s 
can lead to ( 4 Heri) and ( 8 Beri) bound states and can thereby allow for catalyzed 
BBN (CBBN) of 6 Li and 9 Be to abundances far above the ones obtained in standard 
BBN (SBBN) [33, 34, 35, 36]. Indeed, confronting the abundances obtained in CBBN 
with observationally inferred bounds on the primordial abundances of 9 Be (and 6 Li) 
imposes restrictive upper limits of Yf < 2 x 10~ 15 (2 x 10~ 15 - 2 x 10~ 16 ) for > 10 5 s; 
cf. Fig. 5 in Ref. [36] for = n^* . 

For example, in gravitino LSP scenarios with the t\ LOSP being the next-to-lightest 
supersymmetric particle (NLSP) and conserved R-parity, the listed cosmological constraints 
have been confronted with representative values [37, 38] 

Y,* (0.4-1.5) x 10- , (1.1) 

which are in good agreement with the curves in Fig. 1 of Ref. [37] that have been obtained 
for the case of a purely 'right-handed' f\ ~ tr NLSP and a bino-like lightest neutralino, 
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Xi — B, with a mass of = 1.1 m^. Thereby, it has been found that the (C)BBN 
constraints impose the limit r ?1 < 5 x 10 3 s [33, 36] with severe implications in the collider- 
friendly region of < 1 TeV: (i) The limit disfavors the kinematical determination 
of [39] and thereby both the determination of the Planck scale at colliders [16] and 
the method proposed to probe the maximum reheating temperature Tr at colliders [27]. 
(ii) Within the CMSSM, the limit implies an upper limit on the reheating temperature 
of Tr < 10 7 GeV [40, 41, 42] that disfavors the viability of thermal lepto genesis with 
hierarchical heavy Majorana neutrinos [43, 44, 45, 46, 47]. (iii) The limit can point to 
a CMSSM mass spectrum which will be difficult to probe at the Large Hadron Collider 
(LHC) [48, 40, 41, 42]. Indeed, the limit can be relaxed only with a significant reduction 
of (1.1) which has been presented explicitly so far only for non-standard cosmological 
scenarios with a very low value of Tr [49] or with late-time entropy production after t\ 
decoupling and before BBN [40, 50]. 2 

In this work we calculate the decoupling of the lighter stau from the primordial plasma 
by taking into account the complete set of stau annihilation channels in the MSSM with 
real parameters for SUSY spectra for which sparticle coannihilation is negligible. Using 
our own code for the computation of the resulting thermal relic stau abundance Yf, we 
examine explicitly (i) the effect of left-right mixing of the lighter stau, (ii) the effect of 
large stau-Higgs couplings, and (iii) stau annihilation at the resonance of the heavy CP- 
even Higgs boson H°. We consider both the "phenomenological MSSM" (pMSSM) (see, 
e.g., [51]) in which the soft SUSY breaking parameters can be set at the weak scale, and the 
CMSSM, in which the gaugino masses, the scalar masses, and the trilinear scalar couplings 
are assumed to take on the respective universal values rrii/2, mo, and Aq at the scale of 
grand unification Mgut — 2 x 10 16 GeV. Within the framework of the pMSSM, we show 
examples in which Yf can be well below 10 -15 . Even within the CMSSM, we encounter 
regions with exceptionally small values of 1~ < 2 x 10 -15 . The implications of these findings 
are discussed for scenarios with the gravitino LSP and the stau NLSP. We also address the 
viability of a t\-tI asymmetry. Remarkably, we find that key quantities for the significant 
Yf reduction could be probed at both the LHC and the International Linear Collider (ILC). 

A calculation of the thermal relic abundance of long-lived staus has also been part 
of a recent thorough study [52] which focusses on gauge interactions and on the effect of 
Sommerfeld enhancement. In contrast, the most striking findings of our study — in which 
Sommerfeld enhancement is not taken into account — are related to the Higgs sector of the 
MSSM. At this point, we should also stress that the micrOMEGAs code [53, 54, 55, 56] allows 
for sophisticated calculations of the thermal relic stau abundance also in regions in which 
coannihilation effects become important. In fact, micrOMEGAs has already been applied 
in several studies to calculate Yf [38, 40, 57, 42, 52]. In this paper, we also work with 
micrOMEGAs to cross check the results of our own Y? calculation and to calculate Yf in 
parameter regions in which sparticle coannihilations become relevant. 

The outline of this paper is as follows. In the next section we review basic properties 
of the staus to introduce our notations and conventions for the stau mixing angle. Sec- 

2 Note that some implications of the limit can be evaded not only by relaxing it but also by respecting 
it, e.g., R-parity violation can lead to Tf , < 5 x 10 3 s [11]. 
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tion 3 explains the way in which we calculate and provides the complete list of stau 
annihilation channels. In Sect. 4 we analyze the dependence of the most relevant stau 
annihilation channels on the stau mixing angle. Effects of large stau-Higgs couplings and 
stau annihilation at the H° resonance are studied in Sects. 5 and 6, respectively. The 
viability of a t\-t\ asymmetry is addressed in Sect. 7. In Sect. 8 we present exemplary 
parameter scans within the CMSSM that exhibit exceptionally small 1~ values. Potential 
collider phenomenology of the parameter regions associated with those exceptional relic 
abundances and potential implications for gravitino dark matter scenarios are discussed in 
Sects. 9 and 10, respectively. 



2. Stau mixing and mass eigenstates 

In this section we review some basic properties of the stau to set the notation. In ab- 
sence of inter-generational mixing, the stau mass-squared matrix in the basis of the gauge 
eigenstates (tl,tr) reads 

-Ml = H + Tl X \ \ = | "J 1 | R~ T (2.1) 



m T X T m T + m RR 




with 



™LL = m k + + sin2 M l cos 2 P ( 2 - 2 ) 

m | R = m ~ R - sin 2 9 W M% cos 2/3 (2.3) 
X T = A T -/x*tan/3 . (2.4) 

Here, and m^ R are the soft SUSY breaking masses, A T is the trilinear coupling, \i is 
the Higgs-higgsino mass parameter, and tan/3 = V2/v\ denotes the ratio of the two Higgs 
vacuum expectation values. In this work we restrict ourselves to the MSSM with real 
parameters. Then X* = X T so that the mass eigenstates t\ and T2 are related to tl and 
tr by means of an orthogonal transformation 

. , _ / cos#= sinfe \ ,„ „. 

with R 7 =[ . T J 2.5 

\ — sin % cos % / 

with 9^ denoting the stau mixing angle. Imposing the mass ordering < m^ 2 and 
choosing < 8? < it, the mixing angle can be inferred from the elements of 

2m T X T 2m T X T . 2m T X T 
tan 28r = — n 5— = : — , sin 26^ = — 5 5- , (2.6) 

m LL ~ m RR " m $ 1 — m ?2 

where the sign of the second relation determines the quadrant of 8^. In the first relation, 
we have introduced S = m 2 ^ — ?n RR . In particular, 9? = ir/2 corresponds to a purely 
right-handed stau, T\ = tr, whereas maximal mixing occurs for 8^ = 7r/4 and 37r/4. The 
physical stau masses are then given by 

1 

2 




m n 2 = m r + m RR + 



5 T \A 2 + ^m 2 T xf\ (2.7) 



from which we see that an increase of \X T \ leads to a reduction of m^. 
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3. Calculation of the thermal relic stau abundance 



We have undertaken the effort to set up our own full-ffedged relic abundance calculation. 
Let us in the following give a description of our approach to compute the stau yield Yf. 
Throughout this work we assume a standard cosmological history with a temperature T of 
the primordial plasma above the stau decoupling temperature Tf so that the lighter stau 
T\ was once in thermal equilibrium. Then, the total stau yield 1~ = l~ x + Yj» is found by 
solving the well-known Boltzmann equation 

^ = s(crv) [y~ 2 - (Y? q ) 2 ] . (3.1) 
Using the Maxwell-Boltzmann approximation, the stau equilibrium yield Y~ q is given by 

Yp = ^K 2 (^) (3.2) 



7T 2 S V T 

and the thermally averaged annihilation cross section by 



<<™>( T ) = ,^ ds ^S Kl (^j P^a(s) , (3.3) 



where Ki is the modified Bessel function of order % and P e ff = y s — 4m~ ^ 2. 

Note that (<td) contains all the information from the particle physics side. It is obtained 
by computing the total stau-annihilation cross section, 

<r = ^ °tot with cr to t = <*n n-rr + X] a?r i T? ^ ' ( 3 - 4 ) 

where the sum for the annihilation of t\ t± pairs 3 has to be taken over all final states X. The 
factor 1/2 is convention but gives (3.1) its familiar form. The complete list of annihilation 
processes in the MSSM with real parameters — save for coannihilation processes — is given in 
Table l. 4 In addition, this table shows all possible particle exchanges, where s, t, and u are 
the Mandelstam variables which denote the respective channel. A number of annihilation 
processes proceeds also via a four-point vertex. Those are marked in the column named 
"contact." Already by mere optical inspection, we immediately see that the Higgs sector 
plays potentially an important role in the determination of the stau yield Yf . 

For all channels in Table 1, we generate Fortran code for the squared matrix elements 
| Mi\ 2 by using the computer algebra packages FeynArts 5 . 4 [60, 61] and FormCalc 5 . 3 [62, 
63]. For a chosen point in the SUSY parameter space, we then compute the radiatively 
corrected superparticle spectrum by running the spectrum generator SuSpect 2.40 [51]. 
Its output allows us to set all SUSY parameters so that we can compute the total cross 



3 Counting wise we distinguish between f\ t* — * X and the conjugate process t* ti — > X. In absence 
of CP violation in the SUSY sector, their cross sections agree so that we can solve a single Boltzmann 
equation (3.1) for obtaining Yy. 

4 For a purely right-handed stau t\ — tr, the stau annihilation channels and associated cross sections 
have already been presented in Ref. [59] in the context of x?-ri coannihilation. 
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Table 1: The complete set of stau annihilation channels in the MSSM with real parameters for 
scenarios in which sparticle coannihilations are negligible. The mass eigenstates of the Higgs fields 
are denoted by h°, H a , A , and H ± and the ones of the neutralinos, the charginos, and the tau 
sneutrino by Xi .. 4; XT~2! an d v T , respectively. Because of the absence of a tiTiA coupling (cf. 
Sect. 5), s-channel exchange of the CP-odd Higgs boson A and also T\f* — > jA° do not appear. 



(*) (*) 


final state 


s-channel 


i(w)-channel 


contact 




TT (ri ) 




~o 

Xl,..,4 




T\T X -> 


final state 


s-channel 


i(it)-channel 


contact 




fip, ee 


h°,H°, 7 ,Z 


— 


— 




TT 




Xl,..,4 


— 




v e V ei Vfjpn 


Z 




— 




v T V T 


z 




— 




Hk<4k 














Tl 


/ 




zz 






/ 




w+w- 






/ 




7 /i°, 7 i7° 


— 


n 


— 






Z 


Tl, 2 






ZA° 


h°,H° 










h°,H° 








h°h°, h°H°, 
H°H° 


h°,H° 


Tl, 2 


/ 




A°A° 


h°,H° 


T2 


/ 




h°A°, H°A° 


Z 


T2 






H+H- 


h°, H°, 7, Z 




/ 


t k = u, d, c, 


s, t, b 









section Otot(s) given by (3.4) and subsequently the thermally averaged cross section (3.3). 
Numerically, the computation of (3.3) is the most demanding part in the relic abundance 
calculation. In particular, we take special care about the following cases: 

• i?°-resonance: Resonant stau annihilation via H° exchange is one of the central 
points in this paper. In the generation of the matrix elements, we have therefore 
included the total i7°-width T^o in the respective s-channel propagators. 

• Propagator poles: A diverging t(u)-channel propagator can be encountered when a 



- 6 - 



production threshold is met. We overcome this problem by including a "sparticle- 
width" of O.Olm^ in the respective propagators in the vicinity of dangerous thresh- 
olds. A particularly interesting example with a diverging t(it)-channel propagator is 
given by the process t\ r* — ► jH° if y/s = m H o is fulfilled since then the if°-exchange 
in the s-channels of other processes is resonant simultaneously. 

• Bessel functions: The Bessel functions in (3.2) and (3.3) exhibit an exponential be- 
havior for large arguments i»l [64] 



For small temperatures T, the arguments of K\ and K<i in (3.3) become large simul- 
taneously. Therefore, in order to ensure numerical stability, we expand the Bessel 
functions in (3.3) for m^/T > 35 as in (3.5) and cancel the exponents analytically. 5 

We find the starting point for the numerical integration of (3.1) by solving [54] 



where g*(T) is a degrees of freedom parameter [58] and Mp = 2.4 x 10 18 GeV the (reduced) 
Planck mass. Tfi marks the point at which the stau starts to decouple chemically from 
the background plasma, l~(T f i) - F~ q (T f i) ~ \Y~ q (T {1 ) with A = 0.1 [54] chosen in our 
code. Since we use a globally adaptive Gaussian integration routine to calculate (3.3), 
the computation of (crv)(T) is time-demanding. Therefore, we evaluate (3.3) on a grid of 
different temperatures and use cubic spline interpolation to obtain values in between. We 
then solve the Boltzmann equation (3.1) by numerical integration from Tfi to zero. There, 
we fully take into the account the temperature dependence of g* and g*s by interpolating the 
respective tabulated values provided as part of the relic density code DarkSUSY 4.00 [65]. 
The freeze out temperature can then be defined by Tf = (Tfi + T2) /2 where Tf2 is given by 



Y~ q (Tf2) = !%(Tf2)/10 [54]. For T < residual annihilations will further reduce Yf so 



For simplicity, we call this yield Y? in the Introduction and in the following. Moreover, we 
will quantify T in terms of x = m^/T and in particular Tf in terms of Xf = m^/Tf. 

Note that we have additionally modified the FeynArts MSSM model file for the gen- 
eration of the matrix elements in two ways: The first version, which we use throughout 
Sects. 4-6, allows us to set all c/fc^-Higgs and all trilinear Higgs couplings by using the 
computer tool FeynHiggs 2.6.3 [66]; see also Sects. 5 and 6. The second version allows 
for a direct comparison with the existing computer code micrOMEGAs 2.0.6 [53, 54, 67]. 
We have transcribed their routine [68] for the computation of the running quark masses 
to Fortran, adopted all qkQk "Higgs couplings, and modified all Higgs-self couplings of our 

5 J. P. is grateful to P. Gondolo and J. Edsjo for pointing out that trick. 




(3.5) 




(3.6) 



that we refer to the decoupling yield Y~ 



as the quantity at the endpoint of integration. 



-7- 



matrix elements to match with their implemented version of the MSSM [69]. Using this 
second version, we find perfect agreement between our codes. 6 

4. Dependence of stau annihilation on the stau mixing angle 

In order to isolate the distinct features of the different annihilation processes we need 
to have full control over the superparticle mass spectrum. Therefore, in the following, 
we will not rely on any constrained model (such as the CMSSM) where the soft-SUSY 
breaking parameters are subject to stringent boundary conditions at some high scale (such 
as Mgut)- In those models, the mass spectrum is found only after renormalization group 
(RG) evolution from the high scale down to the electroweak scale. Instead, we choose 
to work in the framework of the "phenomenological MSSM" (pMSSM), see, e.g., [51]. 
There, all soft-SUSY breaking parameters can be set at the scale of electroweak symmetry 
breaking — a low scale — which we fix to ~ 1m^ x . In particular, one can also trade the Higgs 
mass-squared parameters m 2 Hu and rn 2 Hd against [i and the pseudoscalar Higgs boson mass 
m A o. 7 Choosing ji as an input parameter is very convenient for two reasons: First, together 
with the specification of the gaugino masses we have control over the gaugino/higgsino 
mixture of the neutralinos Second, n enters directly into the stau-Higgs couplings, 
whose importance will become clear in the next section. Furthermore, in the following, 
we choose to set all soft-SUSY breaking scalar masses (apart from m^ L and m^ R ) to a 
common value Ms = 1 TeV. Thereby, we essentially decouple all sfermions which are not 
of interest for us. This ensures also that we never enter accidentally any coannihilation 
regime. Finally, for simplicity, we set also all trilinear parameters to a common value A. 
Given fj,, A T = A, and tan/3, and thereby X T , we can then fix and 6^ to arbitrary 
values by adjusting m^ R and 5 in Eqs. (2.6) and (2.7). 

In the following, we will focus on two distinct regions of the SUSY parameter space. In 
the beginning, we will choose m A o to be very large m A o = 1 TeV S> Mz. This corresponds 
to the decoupling limit of the MSSM where the following (tree-level) relations hold [70] 

m 2 h0 ~ Af| cos 2 2/3, 
jtih± = m A o + Myj, 

Therefore, m A o ~ m H o ~ m H ± up to corrections O {M\jm A o) so that any of the stau 
annihilation channels into heavy Higgs bosons is kinematically blocked. Furthermore, 
cos (/3 — a) = up to corrections O {M^/m 2 A0 ) implies that the H°VV coupling (V = 
Z, W) becomes very small so that we loose the i7°-exchanges in the stau annihilation 

6 For our computation we use the Standard Model parameters m t = 172.5 GeV, m b (m b ) MS = 4.25 GeV, 
af^iMz) = 0.1172, Q e "™(M z ) = 127.932, and M z = 91.187 GeV. Since micrOMEGAs has hard-coded 
sin 9w = 0.481 from which it computes Mw using the on-shell relation with Mz, we follow their convention 
to allow for a better comparison of our results with micrOMEGAs. 

7 Though the advocated procedure may require fine-tuning in the electroweak symmetry breaking con- 
ditions, it conveniently provides us with running parameters at the scale of stau annihilation. 



m 



i 2 H o ^ m 2 A0 + Mi sin 2 2/3, 

2 M| sin 2 4/3 
cos (/3-a) ~ n 4 . 



(4.1) 
(4.2) 
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Figure 1: (a) The dependence of m? 2 (curved solid line) and m^ T (dashed line) on the stau mixing 
angle 0? for the input parameters = 130 GeV (horizontal solid line), tan/3 = 10, m^a = Ms = 
M 3 = -A = 1 TeV, and 6Mi = M 2 = (J, = 1 TeV (i.e., ^ 5) for which m^o = 169 GeV (dash- 
dotted line) and m^a =116 GeV (dotted line), (b) Dominant stau annihilation cross sections times 
the relative velocity u r of the incoming staus as a function of Or for P e g = 10 GeV and the same 
input parameters as in (a). The curves show the channels with the following final states: h°h°, 77, 
tt, WW, tt, ZZ, -fZ (at Or — 40°, from top to bottom). In addition, we plot a TT Vr for the case of 
a wino-like neutralino, Xi — W, with m^a = 175 GeV as obtained with Mi = 6M2 = 1 TeV (thin 
gray line). No lines are shown for 0^ < 18° where < m^. 



channels with a VV final state. At the same time, the light Higgs boson takes on its Stan- 
dard Model value for the hPVV coupling. Complementary to that we will consider also 
regions of the SUSY parameter space with smaller m^o, e.g., in the next section, where we 
will put a stronger focus on the Higgs sector and its connection to Yf. 

In Fig. la we show the ^-dependence of the masses of the heavier stau, m^ 2 , (curved 
solid line) and the tau-sneutrino, rriy T , (dashed line) for fixed m= = 130 GeV and the input 
parameters tan (3 = 10, m^o = fx = —A = 1 TeV, and 6M1 = M23 = 1 TeV. Because of 
SU(2) gauge invariance, sets also the soft-breaking mass for the tau-sneutrino hence 
approximately m~ ~ + 5 so that v T becomes lighter than tl for 9^ < 18° (6 is negative 
in that region). In addition, we plot the masses of the lightest neutralino, m^o = 169 GeV 
(dash-dotted line), the lighter stau, m= = 130 GeV (horizontal solid line), and the lightest 
Higgs, m h o = 116 GeV (dotted line). We note in passing that m~ may deviate slightly 
from its anticipated input value due to radiative corrections. We then correct for this by 
an adjustment of m~ R so that we indeed ensure to be constant. 

In Fig. lb we plot the dominant stau annihilation cross sections times the relative (non- 
relativistic) velocity in the center-of-mass frame of the incoming staus, v T = 2P c s/m T \, for 
the same parameters as in Fig. la. Owing to an (approximate) Maxwell-Boltzmann distri- 
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bution of the stau velocity, (P e s)\T f ~ y^m^Tf, we choose P e ff = 10 GeV as a representative 
value. 8 The curves show the annihilation channels with the following final states: h°h°, 
77, rr, WW, tt, ZZ, 7Z (at 0? = 40°, from top to bottom). All channels except 77 
show a strong dependence on 0?. The h°h° (tt) channel peaks at 0? = 7r/4 — a feature 
which we will discuss in detail in Sect. 5. For the rr channel, the overall size of the cross 
section is governed by m~o since this channel proceeds only via t(u)-channel exchanges of 
neutralinos. Our chosen input values lead to a bino-like neutralino, x^ 1 — B, and a TT drops 
for an increasingly 'left-handed' stau. (For comparison, the thin gray line shows cr TT v T for 
the case of a wino-like lightest neutralino, = W, of similar mass, m^o = 175 GeV, 
as obtained by changing the gaugino mass input parameters to Mi = 6M2 = 1 TeV.) 
The annihilation into a WW pair becomes important for an increasing tl component in 
ri, i.e., towards smaller 0?, since the t(u)-channel exchange with the tau-sneutrino opens 
up; the T\v r W (j{t\WW) coupling is proportional to cosO? (cos 2 67^). The modulation 
of the 7Z channel can be understood by considering the structure of the ririZ coupling 
oc (1 — 4 sin 2 9 w + cos 20?). Note that the first two terms practically cancel out. For stau 
annihilation into a ZZ pair there is an additional contribution from ^-exchange with the 
respective T1T2Z coupling oc sin 20?. Having discussed the dominant 77 annihilation chan- 
nels in a simple manner, we also warn the reader that interferences between the different 
Feynman diagrams of a given channel may well lead to a counterintuitive behavior. In 
this regard, see Ref. [52] for a thorough discussion of t\t\ annihilation into vector bosons. 
For the limiting case of a purely 'right-handed' stau, t\ ~ tr (0? — > 7r/2), we recover the 
relative importance of the annihilation cross sections into 77, 7Z, ZZ, and rr with bino 
t(it)-channel exchange found in Ref. [37]. 

Figure 2 shows the ^-dependence of Y? (upper panel) and of the relative importance 
of the dominant thermally averaged cross sections, (&iv) / '{(Ttot v ) , at x = 25 (lower panel) 
for the same input parameters as in Fig. 1. The lines in the lower panel are associated with 
the same dominant annihilation channels as in Fig. lb. In addition, the relative importance 
of the sum of the displayed cross sections, {cr^spv) / (a to tv) , (thin line, as labeled) is shown 
to demonstrate that the displayed channels constitute indeed (up to at most about 10%) 
the dominant part of (atotv) for the chosen set of input parameters. In the upper panel, 
the total stau decoupling yield obtained with our own relic abundance calculation is shown 
by the thick line and the one computed with micrOMEGAs, Y~^, by the thin gray line. 
For 0? < 25°, both curves start to deviate from each other since one enters the 17-77 
coannihilation region in which the stau decoupling yield increases. This coannihilation 
effect leads also to the rise of the thin gray line that shows (otot^) Y™ n in arbitrary units 
(a.u.) in the lower panel. Note that the same line illustrates Y? oc l/(<7tot u ) f° r 0? > 25°, 
where the result of our relic abundance calculation agrees with Y- 1 ^. Interestingly, for the 
given input parameters, Y? is not overly affected by the variation in 0? in this region, which 
reflects the fact that ((J tot v) and thereby (av) vary by less than a factor of about 1.5 at 
the relevant time of decoupling. In the next sections, we will demonstrate that this picture 



8 This value is actually at the somewhat lower end, given > 100 GeV and Tt ~ m^/25. However, 
av r depends only weakly on P g, and the thermally averaged (aiv) will be shown in the upcoming figures. 
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Figure 2: Dependence of the stau yield Y? (upper panel) and of the relative importance of the 
dominant thermally averaged cross sections, (<Jiv) / (atotv) , at x — 25 (lower panel) on the stau- 
mixing angle 9^ for the same input parameters as in Fig. 1. In the upper panel, the thick line 
shows the stau yield Y? obtained with our relic abundance calculation and the thin gray line the 
one obtained with micrOMEGAs to which we refer as Y~ an . In the lower panel, the line styles are 
associated with the same dominant annihilation channels as in Fig. lb. In addition, we show (as 
labeled) the relative importance of the sum of the displayed cross sections, (<7di sp f) / (ctotv) , and 
{o'totv} Y- 1 ^ in arbitrary units (a.u.). No lines are shown for 9^ < 18° where rriy T < m^. 

changes significantly for certain other choices of the input parameters. 
5. Effects of large stau-Higgs couplings 

Owing to the scalar nature of the stau, there exists a remarkable difference between the 
standard neutralino decoupling and the scenario in which the long-lived stau freezes out 
from the primordial plasma. For the neutralino LSP, the \i parameter enters into the 
annihilation cross sections only indirectly by influencing the gaugino/higgsino mixture of 
Xi- This stands in strong contrast to the case in which a scalar particle is the lightest 
Standard Model superpartner: the sfermions couple directly to dimensionful parameters 
of the theory, namely, the trilinear couplings A and the Higgs-higgsino mass parameter \i. 
The corresponding operators in the MSSM Lagrangian always contain a Higgs field. In 
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particular, the stau-Higgs couplings are given by 

£MSSM 9 J^— ^2 ^aC[T*,Tp,H]TpH (5.1) 
W a,B=L,R 

with H = h°, H°, A . We have pulled out the factor <?/M\y so that the 'reduced' couplings 
C[t*,T/3,7{] among the gauge eigenstates tl and tr are given by [71] 



C[t*,t, h° 



2 J cp 2 \ cp cpj 

— A T h li— -s w M z s a+f 3 + m T — 

2 V C B CBj CpJ 



(5.2) 



(Tfl \ 
+i— (A T t&n(3 + u)\ 

+ 2 „ J ' < 8S > 

where C[r*, r, H°] can be obtained from (5.2) upon the replacement a —> a — it/2. When- 
ever convenient, we use the shorthand notation Syy = sin 2 9w, c 7 = cos 7, and s 7 = sin 7. 
The parameters A T and \i only appear off-diagonal and they are multiplied with the asso- 
ciated fermion mass, the tau mass m T . 

Using C = R^CEX, one obtains the couplings of the mass eigenstates t\ and r-i. In this 
regard, it is important to note that the coupling of the CP-odd Higgs boson to the lighter 
stau vanishes, C[r*, 77, A ] = 0. Therefore, we have not listed the process t\t* — > 7 A 
in Table 1. By the same token, there is also no s-channel exchange of ^4° in any of the 
annihilation channels. Note that this statement remains valid even after the inclusion of 
radiative corrections: There is no induced mixing between h°(H°) and A in absence of 
CP-violating effects in the SUSY sector. 

Let us now turn to the probably most interesting couplings in the context of t{t\ 
annihilation, namely, the ones of the lighter stau to h° and H°. The 'reduced' TiTih 
coupling reads 

This is a complicated expression. However, if we choose m^o to be large, m^a 3> Mz, we 
can simplify (5.4) by using cos (j3 — a) = [cf. (4.2)], 



(5.5) 



Thereby, we make an interesting observation: In the decoupling limit (DL), the TiT\h° 
coupling becomes proportional to the left-right entry m T X T of the stau mass-squared 
matrix (2.1) and to sie~- Therefore, it comes as no surprise that the t{t\ annihilation 
cross section into h°h peaks at 8^ = 7r/4 — the point of maximal tl-tr mixing — as can 
be seen, e.g., in Fig. lb. Analogously, one finds that the tiT\H coupling is proportional 
to (A T tan (3 + fx) s 2 e~ in the decoupling limit. Complementary, the TiT2h°/H° couplings 
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exhibit in this limit the same combination of A, fi, and tan (3 as their t{t\ counterparts but 
those terms are now multiplied by C2e~ instead of S20~- 

After the above discussion, it is clear that there exists the possibility to enhance the 
total stau annihilation cross section <7 tot — and thereby to decrease Yf oc l/(cr tot ?;) — by 
choosing a proper combination of large A, fi, and tan/?. In the remainder of this section, 
we will explore this possibility for two exemplary pMSSM scenarios. 

Before proceeding let us make some technical comments. Large values of the pre- 
viously mentioned parameters may well lead to large radiative corrections. 9 In order to 
arrive at a proper loop-improved tree-level result, we re-evaluate the entire Higgs sector 
using FeynHiggs. In particular, we have modified our generated matrix elements in a way 
that allows us to set all trilinear Higgs couplings to their loop-corrected values. 10 Note 
that this goes well beyond a simple a — > a e s prescription. Only then, we mostly find 
better agreement of our cross sections for stau annihilation into two Higgses with the ones 
computed by micrOMEGAs. The latter program uses CalcHEP [68] for the generation of the 
matrix elements. There, the trilinear Higgs self-couplings have been expressed in terms of 
m h o, m H o, and m A o which effectively reabsorbs a bulk of the radiative corrections [69]. 
We therefore think that we do slightly better whenever we encounter some disagreement 
between the mentioned cross sections. Though the overall effect on Yf is typically small, 
it can be at the level of 20% (see below). Finally, it is well known that a large A param- 
eter may lead to charge/color breaking minima (CCB) in the scalar MSSM potential; see, 
e.g., Ref. [73]. SuSpect performs some basic checks which we take into account to make 
sure that we do not violate the constraints associated with CCB. In any case, our pMSSM 
scenarios shall be regarded as toy models which allow us to extract important features of 
primordial stau annihilation in the most transparent way. 

In Fig. 3 we demonstrate the effect associated with a large T\T\h° coupling by presenting 
the ^-dependence of Yf (upper panel) and of the relative importance of the dominant 
thermally averaged cross sections, {(Jiv) / {a to tv) , at x = 30 (lower panel) for the pMSSM 
scenario associated with = 130 GeV, tan (3 = 50, m A o = Ms = M3 = —A = 1 TeV, 
and 6M1 = M2 = fi = 1 TeV. In this scenario, m h o stays in the range 117 — 119 GeV 
and the lightest neutralino is bino-like with a mass of m^o = 169 GeV. Stau annihilation 
into heavy Higgses remains kinematically forbidden. The curves in the lower panel are 
associated with stau annihilation into h°h°, WW, tt, ZZ, 77, and 7Z (at 9^ = 80°, 
from top to bottom). As is evident, the annihilation into h°h° is enhanced already well 
before 6^ = 7r/4. At the peak position, <j h a h ov T ~ 8.8 x 10 3 pb for P c g = 10 GeV (no 
thermal average), which is still three orders of magnitude below the unitarity bound for 
inelastic s-wave annihilation, cr n v T = 87r/(m^ 1 P e ff) [74, 52]. Also the cross sections for 
stau annihilation into WW and ZZ are strongly enhanced towards % = 7r/4 since the 
s-channel contribution of t{t\ — > h°* — ► VV becomes very important. At their respective 
peak positions, awwVr — 250 pb and azzv T ~ 130 pb for P c r = 10 GeV. (Because of the 
dominance of the h°h° channel, the corresponding maxima do not show up in Fig. 3 where 

9 In this context, note that we introduce a large mt-m^ splitting when choosing Ms = 1 TeV. 
We are grateful to T. Plehn and M. Rauch for providing us, for cross-checking, with their implementation 
of a Fortran routine which calculates the Higgs self-couplings using the effective potential approach [72] . 
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Figure 3: Analogous to Fig. 2 but for the pMSSM scenario associated with = 130 GeV, 
tan/? = 50, and m A o = fi = M s = 6Mi = M 2 , 3 = -A = 1 TeV and for x = 30. The stau 
decoupling yield takes on its minimum value of Y? = 7.4 x 10 -16 at 6? = 45°. The displayed stau 
annihilation channels are associated with the following final states: h°h°, WW, tt, ZZ, 77, and 
'jZ (at 6r — 80° , from top to bottom) . No lines are shown for 6? < 4° where m^ T < mf x . 



((Tju) / {(Tt tv) is shown.) By the same token, the cross sections of all (kinematically allowed) 
channels with a fermion-antifermion final state (e.g. tt) — which are subdominant in the 
scenario considered in Fig. 3 — experience an enhancement for 9^ — > tt/4. In total, there is 
an enhancement of (otot^) that delays the thermal freeze out of the staus significantly, i.e., 
Xf ~ 33 for #jr ~ 7r/4. As can be seen in the upper panel of Fig. 3, the decoupling yield is 
thereby reduced dramatically down to a minimum value of If = 7.4 x 10~ 16 for maximal 
left-right mixing of the staus. 

In the previous pMSSM examples, annihilation into final states containing heavy Higgs 
bosons is kinematically forbidden. We can allow for those channels by reducing the input 
value m^o. Indeed, scenarios in which all Higgs bosons are very light in conjunction with 
large tan/3 have been studied in the literature, see, e.g., [75, 76] and references therein. 
We thus consider now the following pMSSM scenario: m^o = 130 GeV, = 150 GeV, 
tan/3 = 50, M s = M 3 = -A = 1 TeV, and 3M 1 = M 2 = fj, = 1 TeV. In Fig. 4, 
the associated ^-dependence of 1~ and of / (cr to tt?) at x = 30 for the now dominant 
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Figure 4: Analogous to Fig. 3 but for the pMSSM scenario associated with = 150 GeV, 
tan (3 = 50, m A a = 130 GeV, M s = M 3 = -A = 1 TeV, 3Mi = M 2 = n = 1 TeV and for x = 30. 
The stau decoupling yield reaches its minimum value of Yf = 4.1 x 10~ 16 at 6? = tt/4. The displayed 
stau annihilation channels are associated with the following final states: h°H°, h°h°, H°H a , bb, 
WW, ZZ, tt, and 77 (at 6? = 75°, from top to bottom). For an optimized presentation of those 
channels, the line indicating the relative importance of the sum of the displayed cross sections is 
scaled down by a factor of 1/2: (odispv) /2(a to tv). No lines are shown for 9^ < 4° where mn T < . 



channels is shown in a similar way as in Fig. 3; only the relative importance of the sum of 
the displayed cross sections is scaled down by a factor of 1/2, (cdisp' u )/2(o"tot^)i to allow 
for an optimized presentation of the single dominant channels. Throughout the considered 
ftp range, the masses of both CP-even Higgs bosons are relatively light and remain rather 
constant: m h a = (118 ± 1.5) GeV and m#o = (128.5 ± 1) GeV. Here the dominant 
annihilation channels are associated with the following final states: h°H°, h°h°, H°H°, 
bb, WW, ZZ, tt, and 77 (at ftp = 75°, from top to bottom). As can be seen, stau 
annihilation into h°H° is now more dominant than the one into h°h° and also the H°H° 
channel becomes important, where each of those channels is indeed associated with an 
(absolute) annihilation cross section {&iv) that peaks at 6^ = 7r/4. Also the annihilation 
into bb is significant — a process which we will discuss in detail in the following section. 
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In this respect, one should stress that all processes with s-channel H exchange are here 
less suppressed by m 2 H0 in the respective propagator than in the previously considered 
scenarios. Note that the asymmetry of {(Tiv) / {a tot v) of those dominant channels (h°H°, 
h°h°, H°H°, bb) with respect to a reflection at 6? = 7r/4 is dominantly caused by the 
^-dependent modulation of the WW channel. As in the pMSSM scenario considered in 
Fig. 3, there is again an significant enhancement of (otot^) that delays the stau freeze out 
such that X{ ~ 33 at 9^ ~ 7r/4. Thereby, the efficient annihilation into final state Higgses 
is accompanied by a significant drop in Y~ down to 1~ = 4.1 x 10~ 16 at 0? = ir/A as can 
be seen in Fig. 4. At this minimum, there is a 20% disagreement between Y? from our 
calculation of stau decoupling (solid line) and the micrOMEGAs result Y- 11 ^ (thin gray line) 
which is a consequence of the different treatments of the Higgs sector described above. 

Let us finally remark that the Higgs couplings to fermions and vector bosons as well as 
the Higgs self-couplings develop a strong dependence on m A o once we leave the decoupling 
regime (m A o < 200 GeV); for a comprehensive review see, e.g., Rcf. [77]. 11 Changes in 
m A o can therefore be accompanied by shifts in the relative importance of the corresponding 
annihilation cross sections. This underlines the fact that the details in the Higgs sector 
may very well be crucial for the determination of the relic abundance of a long-lived t\. 

6. Resonant stau annihilation 

By inspection of Table 1 it becomes clear that primordial stau annihilation can also proceed 
resonantly via s-channel exchange of the heavy CP-even Higgs boson H° for m H o ~ 2m^ . 
While the LEP bound on the stau mass > 82 GeV [1] forbids h° to become on- 
shell (m™ ax ~ 140 GeV, e.g., [77]), the s-channel exchange of ^4° is absent 12 because of 
C[ri,r*, A ] = (see Sect. 5). Again, our choice to work in the framework of the pMSSM 
proves to be very helpful. Since the H° resonance occurs for 2m^ ~ m#o , one runs quickly 
into the decoupling limit in which m H o is governed by the input parameter m A o according 
to the simple relation (4.1). This allows us to scan through the resonance easily. 

Let us explore resonant stau annihilation by considering the exemplary pMSSM sce- 
nario associated with = 200 GeV, 0? = 83° (i.e., a mostly 'right-handed' ti), tan/3 = 
40, and —A = fi = 4Mi = M<i$ = Ms = 1 TeV, for which we vary m A o (and thereby 
m H o) to scan through the resonance. Figure 5 shows the resulting m H o -dependence of Y? 
(upper panel) and of ((Jiv) / {(Ttotv) at x = 25 for the dominant annihilation channels (lower 
panel). Those channels are now associated with the following final states: bb, tt, tt, 77, 
h°h°, WW, >yZ, ZZ, and tt (at m H o = 350 GeV, from top to bottom). In Table 1 all res- 
onant channels can be identified. Close to the resonance condition 2m^ ~ m H o, the most 
important processes are stau annihilation into bb and tt. This is because the couplings 
of those final state fermions to H° are tan/3 enhanced: for tan/3 3> 1, the ffH° coupling 

The Higgs sector is also particularly sensitive to the mixing in the stop sector. In the considered pMSSM 
scenarios, \X t \ = \ A t — /icot/3| ~ Ms which corresponds to the "typical-mixing scenario" [78]. 

12 Even in absence of SUSY- induced CP violation, resonant annihilation via A°-exchange may still proceed 
through T1-T2 coannihilation. However, this scenario requires considerable fine-tuning in the stau mass- 
squared matrix since r\ and T2 have to be nearly degenerate. 



- 16 - 




m H o 



Figure 5: Dependence of Yf (upper panel) and of (aiv) / (<7totv) at x = 25 (lower panel) on m#o 
for the pMSSM scenario associated with m?i = 200 GeV, 0? — 83°, tan/3 = 40, and —A = /i = 
4Mi = Af2,3 = = 1 TeV. In the upper panel, the dark line shows the stau yield Y? obtained 
with our relic abundance calculation and the thin gray line the one obtained with micrOMEGAs. The 
stau decoupling yield takes on its minimum value of Yj= — 9.7 x 10 -16 at m H o = 404 GeV. In the 
lower panel, the displayed dominant stau annihilation channels are associated with the following 
final states: bb, tt, tt, 77, h°h°, WW, jZ, ZZ, and it (at m H <> — 350 GeV, from top to bottom). 



~ mfSp_ a tan (3 with / = b,r [71]. The (broad) peak associated with the resonance 13 
already builds up for m H o > 2m^ = 400 GeV. At zero relative velocity, this would be a 
region in which the H° resonance cannot occur. However, since T\ is in kinetic equilibrium 
at the time of freeze out, resonant annihilation takes place already for 2m^ < mjjo [79]. 
For m H o < 2m^ = 400 GeV, the processes containing s-channel H° exchange proceed 
with a slightly faster rate (if kinematically allowed). The impact of the H° resonance on 
the thermal t\ freeze out and the resulting is substantial. Since the total width of H° 
is F H o = (6 — 10) GeV for m H o = (300 — 500) GeV in the considered pMSSM scenario, 
the reduction of 1% extends over a relatively large m H o range. In this regard, note that 
r^o could be substantially larger had we not essentially decoupled all sfermions — except 

13 Notice that we plot (<Jiv) / (a to tv} so that the actual shape of the resonance looks somewhat different. 
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Ti, T2, and z? r — by choosing Ms = 1 TeV. For m^o ~ 404 GeV, i.e., at the dip of the 
resonance, we find Xf ~ 33 and a minimum stau decoupling yield of Yf = 9.7 x 10~ 16 (dark 
line). Thus, despite the (still) moderate value of tan/3 = 40, a significant reduction of 
is encountered. Indeed, can be even further suppressed for a larger value of tan (5. 
Let us remark that an accurate determination of Y^ in the resonance region requires to 
take special care of the bbH° vertex. This coupling is well known to receive substantial 
radiative corrections for sizable values of tan (3. Therefore, we rely again on the computer 
tool FeynHiggs to compute all quark-antiquark-Higgs couplings and the total width T H o. 
Also the micrOMEGAs code takes special care of the bb~H° vertex. We therefore think that 
the difference between the yields shown in the upper panel of Fig. 5 reflects the theoretical 
uncertainty involved in the determination of T^o as well as the bbH° vertex. 

7. On the viability of a Ti-Tj* asymmetry 

Given the strong bounds on the abundance of negatively charged T\ from bound-state effects 
during BBN, i.e., from CBBN of 6 Li and 9 Be, it is natural to ask whether it is possible to 
have an excess of positively charged r*'s over negatively charged ri's. The generation of a 
particle-antiparticle asymmetry requires a departure from thermal equilibrium. Therefore, 
one might think that a t\-t\ asymmetry can be produced at the time of the stau freeze 
out if the (slepton number violating) process T\j\ — > rr occurs at a different rate than its 
conjugate counterpart. Such a situation might indeed occur if we allow for (CP-violating) 
complex values of the parameters A T , fi, and in the SUSY sector. However, the 

staus are still tightly coupled to Standard Model particles so that they remain in kinetic 
equilibrium with the primordial plasma. Therefore, any excess of r* over T\ arising will 
be washed out quickly by the inelastic scattering process t{t <-> tit. 14 Indeed, it is well- 
known [79] that processes of the latter type occur at much larger rates than the rates for 
the mutual annihilation of the decoupling particle species. The same argument given in [79] 
can be adopted to our case. At the time of freeze out, the reaction rates of interest can be 
estimated as 

Tin -»■ tt : n^n^a^-.TT ~ T z m\e~ 2m ^ /T a^^ TT , (7.1) 
r* T - tit : n^n T ^ T ^- ~ T 9 / 2 m% 2 e~ m ^ T a^ T ^- , (7.2) 

—(*) 

since r{ is approximately Boltzmann distributed. For simplicity, we have treated the tau 
lepton r as a (still) relativistic species. By taking the ratio of (7.2) with respect to (7.1), 

{T/m^f 12 e m? i /T ~ 10 9 for m^/T ~ 25 , (7.3) 

we find that the equilibrating process is by far more dominant. Here, we have used that 

o^ti^tt and OTf* T ^fyf are not too different. In fact, both processes proceed at tree level 

exclusively via Xi exchange so that one cannot decouple (7.2) from (7.1) by a simple 

adjustment of the neutralino mass spectrum. 

14 We are grateful to M. Pospelov for pointing out that process. Additional equilibrating processes are, 
e.g., TiW~ ^> TiW + or t{H~ <-> TiH + , which are however Boltzmann-suppressed. Also note that a 
lepton asymmetry of the order of the baryon asymmetry is expected in baryogenesis scenarios based on 
leptogenesis. 
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Figure 6: Contours of Y? (as labeled) in the (m 1 / 2 ,mo) plane for tan/3 = 43, Ao = 0, and /j, < 0, 
where darker shadings imply smaller Yf values. The dashed lines are contours of = 100, 300, 
and 600 GeV (from left to right). The light-shaded region at rowj < 450 GeV is excluded by the 
LEP bound m^o < 114.4 GeV [1]. In the white area either m^o < or correct electroweak 
symmetry breaking is not established (in the very upper left corner), where the thin contours 
indicate the Higgs funnel in the xl LOSP region. Table 2 provides detailed information for the 
SUSY model represented by the point "A" that is indicated by the star. 



8. Exceptionally small stau abundances within the CMSSM 

We have shown above that the total stau annihilation cross section can be significantly 
enhanced. The thermal freeze out of ri's is thereby delayed such that their abundance 
prior to decay, Yf, is suppressed. In the following we focus on the CMSSM to see whether 
the effects discussed in Sects. 5 and 6 do appear also in models in which the pattern of 
soft-SUSY breaking parameters fulfills certain boundary conditions at a high scale. Note 
that we compute Yf with micrOMEGAs in this section since coannihilation processes are 
not included in our relic density code. In addition, we employ SPheno 2.2.3 [80] for 
the computation of the mass spectrum and the low energy constraints associated with 
B{b — > 57) and the anomalous magnetic moment of the muon Let us now proceed by 
discussing two exemplary CMSSM parameter scans. 

Figure 6 shows contours of constant Yf in the (m 1 /2,mo) plane for tan (3 = 43, A$ = 0, 
and a negative sign of the \i parameter. The contour lines represent the values Yf = 10~ 14 , 
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tan/3 = 55, Aq = 2mo, /x > 
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Figure 7: Contours of Yf (as labeled) in the (mi/ 2 ,mo) plane for tan/3 = 55, Ao ~ 2mo, and 
/i > 0, where darker shadings imply smaller Yf values. The dashed lines are contours of = 
100, 300, and 600 GeV (from left to right). The large light-shaded region in the lower left corner 
is excluded by bounds from direct Higgs and SUSY searches (or by the appearance of a tachyonic 
spectrum). In the region to the left of the vertical solid and dotted lines, m^o < 114.4 GeV [1] 
and B(b — * sj) > 4.84 x 10~ 4 [81], respectively. In the white area, m^o < m^. Table 2 provides 
detailed information for the SUSY models represented by the stars "B" and "C" (as labeled). 



4 x 10 , 10 , and 4 x 10 , where darker shadings imply smaller values of Y?. The 
dashed lines are contours of = 100, 300 and 600 GeV (from left to right). The light- 
shaded region at m^j < 450 GeV is excluded by the mass bound m h o > 114.4 GeV from 
Higgs searches at LEP [1]. The white area indicates the region in which either correct 
electroweak symmetry breaking is not established (in the very upper left corner) or in 
which m~o < m^. Since \i < 0, the plane is actually in tension because of (negative) SUSY 
contributions a^ USY to the anomalous magnetic moment of the muon, a^ = (g — 2)^/2. 

Figure 7 presents a scan over the (m^2, mo) plane for tan (3 = 55, Aq = 2mo, and fi > 
with contours of 1~ = 4 x 10~ 15 , 10~ 14 , 4 x 10 -14 , 10~ 13 , and 4 x 10~ 13 (darker shadings 
indicate smaller Yf values) and = 100, 300, and 600 GeV (dashed lines, from left 
to right). The large light-shaded region in the lower left corner is excluded by the robust 
bound > 82 GeV [1] from collider searches of charged sleptons (or by the appearance of 
a tachyonic spectrum). The LEP Higgs bound m h o < 114 .4 GeV [1] is situated within this 
region in close vicinity to its boundary for mo < 400 GeV and is indicated by the solid line 
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for m > 400 GeV. In the region to the left of the dotted line, B(b -> 57) > 4.84 x 10~ 4 [81], 
which is in tension with the bounds from inclusive b — ► 57 decays. 

Let us now discuss some generic features of the stau yield within the CMSSM on 
the basis of Figs. 6 and 7. We note beforehand that our more general statements on the 
?i LOSP region in the CMSSM are corroborated by a parameter scan over the following 
range 15 

m 1/2 = (0.1 - 6) TeV, tan/3 = 2- 60, 

- 4m < A) < 4m , sgn// = ±l. (8.1) 

In both figures an almost horizontal, narrow band of low Y~ appears in which 2m^ 1 ~ m^o 
holds so that stau annihilation proceeds via resonant production of the heavy CP-even 
Higgs boson H°. We have marked points the centers of the respective regions with "A" 
and "B" for which we provide detailed information in Table 2. Given a present uncertainty 
of ~ 3 GeV in the determination of m h o [82], we note that the LEP Higgs bound has to 
be treated with some care. For example, a (vertical) m h o = 112 GeV contour would be 
situated at m 1 / 2 — 400 GeV in the resonance region of Fig. 7. Accordingly, one could 
consider the entire resonance region shown to be compatible with direct Higgs searches. 
However, due to the large value of tan (3 = 55, the bound on b — ► is very strong so that 
a large part of the resonance region remains excluded by this constraint. In this regard, 
it is interesting to see (Fig. 6) that 2m^ 1 ~ m#o also appears in the t\ LOSP region for 
lower values of tan (3. In the center of both resonance regions, the yield becomes as low as 
1~ = 4.2 x 10 -15 (point A) and i~ = 2.5 x 10~ 15 (point B). Despite the heavier mass of the 
lighter stau (see Table 2), the suppression of Yf is still more pronounced in Fig. 7 than in 
Fig. 6. This is because the bottom Yukawa coupling becomes larger with increasing tan f3, 
as discussed already in Sect. 6. In fact, annihilation into bb final states is in both cases 
by far the dominant process with relative importances of 76% (point A) and 87% (point 
B). The extension of both resonance regions is due to the total width of H° of respectively 
T H o ~ 9.6 GeV (point A) and T H a ~ 22 GeV (point B); note the logarithmic scales in 
Figs 6 and 7. We note in passing that the appearance of the H° resonance does not imply 
the absence of the neutralino funnel region which is indicated by the (unshaded) contour 
lines in the white area of Fig. 6 

Of course, the question arises whether the appearance of the resonance region is en- 
countered more generically within the framework of the CMSSM. In principle, it is not 
easy to provide a simple quantitative connection between and m H o for arbitrary val- 
ues of the CMSSM parameters. However, without emphasis on an overall applicability, a 
qualitative picture can be drawn. Let us start with the mass of the CP-odd Higgs boson 
m A o which can be written as [83, 84] 

m 2 A0 ~ 1/ sin 2 (3 (ml + 0.52m 2 1/2 + /j? -...). (8.2) 

Here, the ellipsis stand for contributions from the bottom and tau Yukawa couplings. For 
tan/3 > 20, m A0 ~ m\ + 2.5m 2 , 2 — and the corrections from the bottom and tau 

15 Here, we disregard CMSSM parameter points in which SPheno flags an error in the spectrum calculation. 
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Yukawa couplings become important so that m A o is driven towards lower values; note 
that sin 2 /? ~ 1 for tan/3 > 20. Indeed, this property can be used to constrain tan/3 from 
above by confronting m A o with the lower bound from LEP, m A o > 93.4 GeV [1]. On 
the other hand, for large m 1 / 2 , one also enters the decoupling limit of the MSSM so that 
m A o and m H o will be nearly degenerate in mass; cf. (4.1). This can be also seen from the 
exemplary points presented in Table 2. Therefore, also rrijjo will be driven towards lower 
values for growing tan /?. Now, left-right mixing of the lighter stau for not too large values of 
tan /3 is small within the CMSSM, t\ ~ tr, so that approximately m~ ~ mQ+0.15m 2 , 2 [84]. 
Therefore, 2m^ < m H o is the relation that holds usually in the region in which f\ is the 
lightest Standard Model superpartner. However, for large tan/3, the contributions from 
the bottom Yukawa coupling in (8.2) can become strong enough (growing with tuq [84]) 
to overcome any additional decrease of due to left-right mixing so that the resonance 
condition 2m^ ~ m H o can indeed be met. Nevertheless, from scanning over the CMSSM 
parameter range (8.1) it seems to us that the resonance condition 2m^ ~ m H o is not easily 
realized in the part of the t\ LOSP region in which T\-x{ coannihilations are negligible. 
Conversely, it is clear that relaxing the universality conditions for the soft-SUSY breaking 
masses at Mqut will make it easier to find parameter regions in which the resonance 
condition 2ra^ 1 ~ m H o is satisfied. Of particular interest in this respect is the model with 
non- universal Higgs masses (NUHM) with mn 1 ^ rriH 2 m o & t Mqut- There, one can 
adjust the input parameters in order to realize resonant stau annihilation. Indeed, this 
model is qualitatively the same as the class of pMSSM scenarios considered in the previous 
sections, where mn 1 and m# 2 are traded (at the low-scale) against m A o and fj, by using 
the electroweak symmetry breaking conditions. 

Low Yf values are also realized in the narrow vertical region around mi/ 2 ~ 1.1 TeV in 
Fig. 7. At the representative point "C" of that region, Yf = 2.2 x 10~ 15 and the main stau 
annihilation channels are the ones into h°h° (90%) and WW (6%); see Table 2. For larger 
values of m 1 / 2 , Y? exhibits its well known behavior and grows with m^. To the left of the 
Y? = 4 x 10~ 15 contour, the yield increases quickly since the annihilation into h°h° becomes 
kinematically forbidden. Indeed, regions of low Yf which are due to the aforementioned 
annihilation channels are a commonplace appearance in the CMSSM parameter space. 
They are found slightly above the lowest feasible values of m 1( / 2 , i.e., close to the boundary 
of the region which is excluded by direct Higgs and SUSY searches and where > m h o 
still holds. This is because t\ is light in that region since the SUSY particle spectrum 
scales with m 1 / 2 (typically, tuq <C m 1 / 2 for t\ LOSP). Moreover, we find that the LEP 
Higgs bound drops hardly below m 1 / 2 — 450 GeV for tan/3 > 40 and mo < 100 GeV. 17 
Due to a strong correlation between the gaugino mass parameter m\/ 2 and the size of 
the jU parameter, fi 2 ~ (1 — 3)rriy 2 [85], the value of /jl in the experimentally allowed 
region is large. Recall from Sect. 5 that the T\T\h Q coupling is ~ sin26*^V T (m A o » Mz) 
so that \X T \ = \A T — |Utan/3| will become sizeable by increasing tan/?. This leads then 

16 The latter relation ignores contributions from A-terms which can be important but complicate the 
envisaged illustrative picture; for the derivation, we have used m t (m t ) = 163 GeV in Eq. (2.25a) of Ref. [84]. 

17 The position of the LEP Higgs bound (which appears as a near to vertical line for low mo) is very 
sensitive to the value of m t . Lowering m t shifts the bound towards larger values of mi/ 2 . 
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to efficient stau annihilation into h°h° final states. Indeed, in those CMSSM regions, 
also | sin 26^\ is maximized so that Yf already starts to drop below the estimate (1.1) for 
tan/3 > 40. Note, however, that the left-right mixing of t\ within the CMSSM is somewhat 
constrained. Neglecting r- Yukawa contributions, the RG-evolution induced splitting reads 
— ~ 0.37m^ 2 [84] and indeed t\ remains mainly right-handed: By scanning over 
the parameter space, we typically find 65° < % < 115° and thus | sin2#^| < 0.75 in the T\ 
LOSP region in which > m h o and m h o > 114.4 GeV holds. 

9. Prospects for collider phenomenology 

If a SUSY model with a long-lived t\ of < 0.7 TeV is realized in nature, the t\ 
discovery potential will be promising at the LHC with a luminosity of 100 ftr 1 [86]. For 
< 0.25 TeV (0.5 TeV), Ti's can also be examined in precision studies at the ILC 
with a cm. range up to y/s = 0.5 TeV (1 TeV). Once long-lived ti's are produced, one 
should be able to distinguish them from muons by considering the associated highly ionizing 
tracks and with time-of-flight measurements. One should then also be able to infer 
from measurements of the T\ velocity and its momentum [15] and complementary from 
(threshold) studies of the process e + e~ — ► TiT-j 11 at the ILC. 

Remarkably, both mechanisms leading to exceptionally small Yf values come with 
testable predictions: certain ranges of the stau-mixing angle 9^ together with large values 
of tan/3, and/or \A T \ and, in the case of resonant stau annihilation, also m H o ~ 
In particular, the large stau-Higgs couplings lead to an enhanced production of light Higgs 
bosons in association with staus via e + e~ — ► TiT^h and 77 — > rir^/i . The associated cross 
sections can then be relatively large at the ILC with a sufficiently high cm. energy [87]. In 
addition, the above reactions with H° instead of h° in the final state can have also relatively 
large cross sections if H° and t\ are sufficiently light. These reactions will then allow for 
an experimental determination of the stau-Higgs couplings and clarify whether its values 
are compatible with an extremely small value of 1~ [87]. Moreover, a measurement of m H o 
pointing to m H o ~ 2m^ could be an experimental hint for resonant stau annihilation in 
the early Universe. 

Indeed, the scenarios considered could allow for a determination of both m h o and 
m H o already at the LHC. Because of the large values of tan (3, the dominant production 
mechanism for h°/H° will be the associated production of the neutral Higgs bosons with 
bottom quark pairs, pp — > bhh°/H°; see, e.g., [88, 89, 90, 91] and references therein. In 
fact, associated bbh°/H° production with h°/H° — > [i + is considered as one of the most 
promising processes for measurements of m^o at the LHC despite the relatively small 
h° /H° — > ^ + // _ branching ratio [92]. In SUSY scenarios with a sufficiently light long-lived 
Ti LOSP, these processes will be complemented by associated bbh°/H° production with 
h°/H° — > tit*, where measurements of the invariant mass of the t{t{ pair could potentially 
provide a unique way to infer m h o and m H o at the LHC. In fact, h°/H° — > t\t\ will occur 
most prominently exactly in the regions associated with the exceptional Yf values due to 
the enhanced stau-Higgs couplings. Having outlined these proposals, we leave a dedicated 
study for future work. 
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Table 2 illustrates that the kinematical reach of both the LHC and the ILC could be 
sufficiently large to allow for the studies mentioned above. In none of the given points does 

exceed 200 GeV so that t{t\ pair production would already be possible at the ILC 
with \fs < 0.5 TeV. There, one could also produce tit^/i final states in scenarios A and 
C. Even the condition m H o ~ 2m^ could be probed in both scenarios A and B that allow 
for resonant stau annihilation. 

10. Implications for gravitino dark matter scenarios 

As already stressed in the Introduction, 1~ is subject to stringent cosmological constraints. 
Indeed, to decide on the cosmological viability of a SUSY model, one has to confront the 
associated 1~ values with those constraints. In particular, for gravitino LSP scenarios 
with unbroken R-parity, restrictive cosmological constraints and implications thereof have 
been derived [48, 39, 40, 41, 57, 42, 32, 93, 36] often under the assumption that Y~ can 
be described by (1.1). However, while (1.1) is quite reliable for t\ ~ tr [37, 38, 40, 52], 
we have shown in the previous sections that Yf (for a given m^) can be about than two 
orders of magnitude smaller than (1.1). Thus, in gravitino dark matter scenarios with such 
exceptionally small Yf values, our understanding of the cosmological constraints and the 
associated implications could change significantly. 

To demonstrate this point, let us indicate for which >~ values the existing cosmological 
constraints (in their present status) are respected: 

• For 1~ < 10~ 14 , the upper limit on Y~ imposed by the non-thermal production 
of gravitinos in f\ decays, Q ~ TP < / 0<j m — given explicitly in (22) of Ref. [22] — is 
respected for mg < 500 GeV even if only a small fraction / = 0.01 of dark matter 
is assumed to originate from f\ decays; cf. Fig. 13 of Ref. [22]. This applies equally 
to other scenarios with an extremely weakly interacting LSP — such as the axino 
LSP [5, 19] — originating from t\ decays. 

• For 1~ < 10~ 13 , the BBN constraints associated with effects of hadronic energy release 
on the primordial D abundance can be respected for f\ ~ tr and up to 10 TeV 
independent of the t\ lifetime; cf. Fig. 11 of Ref. [22]. For a sizable admixture of tl 
in ti, this Yf constraint can become more restrictive in particular with the enhanced 
stau-Higgs couplings allowing for exceptionally small Yf values. Nevertheless, these 
exceptional values are typically associated with m^ x < 300 GeV where the Yf limit 
is significantly more relaxed: Yf < 10~ n for ti — tr. A tightening to Yf < 10~ 13 
(10 -15 ) will then require an increase of £h a d by a factor of 10 2 (10 4 ). On the other 
hand, sufficiently degenerate mg and will always be associated with small values 
of e had an d thereby with relaxed Yf limits from energy release, even in the case of 
strongly enhanced stau-Higgs couplings. 

• For Yf < 10~ 14 (10~ 15 ), the BBN constraints associated with effects of electromag- 
netic energy release on the primordial D ( 3 He) abundance can be respected indepen- 
dent of the Ti lifetime; cf. upper panels of Fig. 12 (100 GeV < < 10 TeV) of 
Ref. [22] and Figs. 14 (m 7l = 100 GeV) and 15 (m 7l = 300 GeV) of Ref. [32]. 
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• For Y~ < 2x 10~ 15 (2x 10~ 16 -2x 10~ 15 ), the BBN constraints associated with bound 
state effects allowing for CBBN of 9 Be ( 6 Li) can be respected even for > 10 5 s; 
cf . Fig. 5 in Ref . [36] . The uncertainty on the 1~ limit associated with CBBN of 6 Li 
reflects the difficulties in quantifying an upper limit on the primordial 6 Li abundance; 
for more details on this issue, see [36] and references therein. Note that the given 
limits correspond to upper limits on the primordial fractions of 9 Be/H and 6 Li/H of 
2.1 x 1CT 13 and 1CT 11 - 1(T 10 , respectively. 

Thus, the SUSY models which come with thermal relic stau abundances of Yf < 2 x 10~ 15 
can respect each of these cosmological constraints independent of the stau lifetime if a 
primordial 6 Li/H abundance of about 10~ 10 is viable. In particular, the limit of < 
5 x 10 3 s and its implications — discussed in the Introduction — are then no longer valid even 
for a standard cosmological history with primordial temperatures of T > Tf. Thereby, 
the regions with Y~ < 2 x 1CT 15 are associated with particularly attractive gravitino dark 
matter scenarios: 

• The gravitino mass can be within the range 0.1 < mg < for which its kinematical 
determination could be viable [16, 21, 24]. Together with measurements of m^ x and 
Tf 15 a kinematically determined would allow one to measure the Planck scale Mp 
at colliders [16, 21, 24]. Indeed, an agreement of the Mp value determined in collider 
experiments with the one inferred from Newton's constant Gn would support the 
existence of supergravity in nature [16]. 

• For mg sufficiently close to m^, the spin 3/2 character of the gravitino becomes 
relevant so that it could be probed in principle by analyzing the decays T\ — > Gt^ [16]. 

• With 1~ < 2 x 1CT 15 , ^g TP is negligible so that basically all of fidm can be provided by 
gravitinos from other sources such as thermal production [25, 27, 28]. Indeed, if also 
gravitino production in decays of scalar fields such as the inflaton [94, 95] is negligible, 
reheating temperatures of Tp < 10 9 GeV could become viable for ~ 100 GeV and 
not too heavy gaugino masses [40, 28]. This would mean that thermally produced 
gravitinos could provide the right amount of dark matter and that thermal lepto- 
genesis [43] with hierarchical heavy Majorana neutrinos — which typically requires 
Tp > 10 9 GeV [44, 45, 46, 47] — would be a viable explanation of the cosmic baryon 
asymmetry, i.e., there would be no gravitino problem. 

• With a kinematically determined iuq, one would be able to probe the reheating 
temperature Tr at colliders and thereby the viability of thermal leptogenesis [27] . 

• For > 10 4 s, the small Yj values could still allow for the primordial catalysis of 
6 Li and 9 Be in agreement with existing astrophysical observations [33, 34, 35, 36]. 

Table 2 illustrates that gravitino dark matter scenarios of the type discussed above 
can even be accommodated within the CMSSM. 18 For gravitino masses of 50 GeV and 

18 Note that this finding points to a caveat of our earlier study [41] relying on < 5 x 10 3 s derived 
from the assumption that is described by (1.1). Our previous sections however show that islands exist in 
which Yr can be significantly below (1.1) even within the CMSSM and for a standard cosmological history. 
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100 GeV, we list the associated values of , of r2~ TP /i 2 , and of the maximum reheating 
temperature T^ ax under the assumption that other gravitino sources can be neglected. The 
expression for can be found, e.g., in (8) of Ref. [41] whereas (12) of the same reference 
can be used to infer the Tg 1 ^ values imposed by ft^h 2 < 0.126. At each CMSSM point 
and for both txIq values, t\ is very long lived, > 10 6 s, and gravitino production from 
T\ decays is negligible, S7~ TP /i 2 < 10~ 4 . In all cases, the gravitino mass mg = 100 GeV 
is sufficiently close to so that the spin 3/2 character of the gravitino can in principle 
be probed [16]. A reheating temperature of Tr > 10 9 GeV is viable only for the points 
A and B with m 1 / 2 significantly below 1 TeV, i.e., at the points at which resonant stau 
annihilation leads to the reduction of Y~. Because of > 10 6 s, the limit from CBBN 
of 9 Be is at V^ ax ~ 2 x 10~ 15 for each point as can be inferred from Fig. 5 of Ref. [36]. 
This bound disfavors point A while the points B and C are associated with values very 
close to this limit and thereby with 9 Be/H ( 6 Li/H) values of about 2.1 x 10~ 13 (10^ 10 ) [36]. 

11. Conclusions 

Supersymmetric models with a long-lived stau T\ being the lightest Standard Model super- 
partner are well-motivated and very attractive in light of potentially striking signatures at 
colliders. For a standard thermal history with primordial temperatures T > m^/20 > Tf — 
which is the working hypothesis is this work — the long-lived t\ becomes an electrically 
charged thermal relic whose abundance can be restricted by cosmological constraints. 

We have carried out a thorough study of primordial stau annihilation and the associated 
thermal freeze out. Taking into account the complete set of stau annihilation channels 
within the MSSM with real parameters for cases with negligible sparticle coannihilation, 
the resulting thermal relic t\ yield 1~ has been examined systematically. While related 
earlier studies focussed mainly on the t\ ~ tr. case [37, 38, 40, 52], we have investigated 
cases in which t\ contains a significant admixture of tl including the maximal mixing case 
and ti ~ tl- 

We find that the variation of the stau mixing angle 6^ does affect the relative impor- 
tance of the different annihilation channels significantly but not necessarily the resulting 
Yf value for relatively small values of tan (3. By increasing tan 0, however, we encounter a 
dramatic change of this picture for large absolute values of the Higgs-higgsino mass param- 
eter /j, and/or of the trilinear coupling A T , which are the dimensionful SUSY parameters 
that govern simultaneously stau left-right mixing and the stau-Higgs couplings: Stau an- 
nihilation into h°h°, h°H°, and H°H° can become very efficient (if kinematically allowed) 
so that Y? can decrease to values well below 10~ 15 . The scalar nature of f\ allows those 
parameters to enter directly into the annihilation cross sections. This mechanism has no 
analogue in calculations of the thermal relic density of the lightest neutralino Xi ■ 

The stau-Higgs couplings are crucial also for the second reduction mechanism iden- 
tified in this work: Even for moderate values of tan/3, we find that staus can annihilate 
very efficiently into a bb pair via s-channel exchange of the heavy CP-even Higgs boson 
H° provided the MSSM spectrum exhibits the resonance condition 2m^ ~ m H o. We have 
shown explicitly that the associated 1~ values can be below 10 -15 as well. This mechanism 
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is similar to the one that leads to the reduction of the Xi density in the Higgs funnel region 
in which neutralino annihilation proceeds at the resonance of the CP-odd Higgs boson A . 

We have worked with an effective low energy version of the MSSM to investigate the 9^- 
dependence of and the two ^-reduction mechanisms in a controlled way. In addition, 
we have shown that the considered effects can be accommodated also with restrictive 
assumptions on the soft-SUSY breaking sector at a high scale. Within the CMSSM, we 
encounter both mechanisms each of which leading to ~ 2 x 10~ 15 in two distinct regions 
of a single (rrii^, mo) plane. 

We have discussed possibilities to probe the viability of the presented l~-reduction 
mechanisms at colliders. While a m H o measurement pointing to m H o ~ 2m^ would sup- 
port resonant primordial stau annihilation, studies of Higgs boson production in association 
with staus, e + e~ (77) — > tit^/i , rirf-ff could allow for an experimental determination of 
the relevant stau-Higgs couplings, for example, at the ILC. Moreover, we have outlined 
that associated bbh°/H° production with h° / H° — > t\t* has the potential to allow for a 
determination of both m h o and m H o at the LHC if a SUSY scenario with large tan (5 and 
large stau-Higgs couplings is realized. 

With the obtained small Y^ values, even the restrictive constraints associated with 
CBBN could be respected so that attractive gravitino dark matter scenarios could be 
revived to be cosmologically viable even for a standard cosmological history. Within this 
class of models, collider evidence for supergravity, for the gravitino being the LSP, and 
for high values of the reheating temperatures of up to 10 9 GeV is conceivable, which could 
thereby accommodate simultaneously the explanation of the cosmic baryon asymmetry 
provided by thermal leptogenesis and the hypothesis of thermally produced gravitinos 
being the dark matter in our Universe. 
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Note added - Ref. [96], in which the potential suppression in the stau yield Y^ due 
to an enhanced annihilation into h°h° final states is also studied, appeared as this work 
was being finalized. This paper provides analytic approximations for the stau annihilation 
cross section into h°h° and for the associated yield. In addition, results of numerical 
studies within the CMSSM, the NUHM, and a scenario with non-universal gaugino masses 
are presented that exhibit parameter regions with extremely small Yf values. In our work 
also enhanced stau annihilation into h°H° and into H°H° and stau annihilation at the 
H° resonance, which were not considered in [96], are discussed. In addition, our work 
provides a systematic investigation based on a complete set of stau annihilation channels, 
an outline of the way in which the mechanisms leading to the suppression of Yjr can be 
probed at collider experiments, and a thorough presentation of the potential implications 
for gravitino dark matter scenarios. 
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Table 2: Exemplary CMSSM points A, B, and C shown in Figs. 6 and 7. In addition to the 
quantities explained in the main text, values of the gluino mass trig and of the mass of the lighter 
stop are given together with the relative importance of the dominant stau annihilation channels, 
x f = m^/Tf, and the decoupling yield Y?. For each point, we list gravitino dark matter scenarios 
with mg = 100 (50) GeV and associated values of the stau lifetime r? 15 the non-thcrmally produced 
gravitino density f2~ TP ft. 2 , and the maximum reheating temperature T^ aK . 
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